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Here, the formalism for the production of fermion-antifermion pairs by vacuum po- 
larisation due to classical fields is investigated according to and extending Refs.^'^. 
The issue of gauge invariance is discussed for the calculation of the expectation 
value for the number of produced pairs from the exact retarded and/or advanced 
propagator. These correlators are given as solutions of the equation of motion for 
the Dirac Green's functions in the classical field. The general characteristics of the 
two point functions in field configurations expected for ultrarelativistic heavy-ion 
collisions are lined out. Calculations for a model for the respective radiation field 
are presented in a covariant-derivative expansion scheme. 



Following the derivation in"^, the expectation value for the number of 
produced pairs can be determined by calculating the overlap of a negative 
energy fermion - corresponding to a positive energy antifermion - with a 
positive energy fermion, after propagation through the field: 

M(p,g) = j d'xd'yi^l{x)GR{x,y)^''^,{y). (1) 

Use has been made of the fact that solutions of the Dirac equation at x and 
y are connected by the retarded propagator which in turn can be expressed 
with the retarded single-particle scattering-operator Tr: 

GR{x,y)^Gl{x~y)+ j d^£,dSGl{x - OrR{^.v)Gl{v - y). 

Free solutions are linked by the free retarded propagator G^j and positive 
and negative energy solutions are orthogonal to each other, whence inte- 
grating the square of the overlap over the phase space (s) yields: 

^/ 2(2t)H 2(2tf3^, ' 

with the unit spinors u{q) and v{p) and where LOp = \/\'p\'^ + ra^ and ujq = 
\/\q\'^ + m^- Gauge invariance can be discussed better in an alternative 
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expression*, where the derivatives act only on the two-point function: 

(n)= J ^^^^dWye^«-(--«)tr{(^/ + m)[z^.-m]G'+_(x,y)[m-i^,]}. 

Even the trace over the correlator G+-{x,y) = {Oin\'il'{y)tp{x) |0i„) is still 
gauge dependent: tv {n{x)G{x,y)n\y)] ^ tr {G{x,y)}. 

I. One way to generate a gauge-invariant result is to project onto gauge- 
dependent solutions in a vacuum configuration of the gauge field [instead 
of free states as in Eq.(l)]. However, the purely partonic interpretation of 
the result is lost in the process: In all gauges other then the one in which 
the gaugc-dcpcndcnt coincides with the free solution, the state contains a 
superposition of an arbitrary number of virtual gluons. Fiirthcr, the inter- 
pretation and the result depend on the choice of that special gauge. This 
method is equivalent to multiplying with adequate gauge (Wilson) links. 
There the choice of the path is the supplementary degree of freedom that 
has to be fixed. Additionally, in the course of the derivation of Eq.(2), q and 
p were interpreted as the momenta of particle and antiparticle, respectively. 
Here, this stays evident only in the one gauge of reference. 

II. If purely fermionic states are to be preserved in all gauges, gauge 
invariance can be achieved by averaging the expectation value over the 
entire gauge group. After a coordinate transformation from d^xd^y to 
d'^[(x + y)/2]d*(x — y) one investigates d{n) /d'^{x — y) and finds that its 
average over the gauge group d{n) /d'^{x — y) is equal to zero but for x = y. 
This is similar to lattice calculations, where the final result contains the 
four-volume V of the lattice as an additional factor. Thus one would obtain 
(n) /V, there. As a; — y is the conjugate coordinate to the fermion momen- 
tum q this implies that the entire dependence on that momentum is only 
due to the geometry of the phase space and the factor ^ -|- m in the trace: 

d{n) _ f '^^1 _Q^^q^-^2T,5{q^ -m'')iT{{i + m)L] 



d^ix-y) J (27r)' 
dH d*k 



L = J ■^^ei-ko)2nSik' -m^)m,m + m)TAik,l) (3) 

; needed retarded and advanced propagators can be written as 
- f d^k ^ ^ 

iGn{x,y,k)^° = ±j -^e'''<^-^e[±{xo - yo)]GH{x,y,k) 
W the homogeneous solution 

GH{x,y,k) = 'Pexp^ J d^ol°[h'dj{x) + ^kj +4.(^0, x) - m]^ 
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Now consider the special case of ultrarelativistic heavy-ion collisions. 
They are frequently described by means of colour charges moving along 
the light-cone^. This picture provides the adequate framework for the de- 
scription of particles with low longitudinal momentum*'^. The classical 
Yang-Mills equations have to be solved for that current. In Lorenz gauge 
the field consists of Weizsacker- Williams sheets and the radiation field in- 
side the forward light-cone. The former represent the Coulomb fields of 
the charges boosted into the closure of the Lorentz group, whence they are 
(5-distributions in the light-cone coordinates: 

Nl,r 

with 27ra^^^^ = ^9ia{^a'^)nL,R lii-^I^T — ^)"i.,rI ^"^^ where A regularises 
the logarithm but does not appear in the field tensor, arc the generators 
of S'C/(3)c. {ta'^)nL,R represent the colour of the charges. It should be 
mentioned that the colour neutrality of each nucleus has to be ensured and 
that the charges of the first nucleus precess in the field of the second and 
vice versa, causing further shcct-likc> contribTitions that can be incorporated 
by generalising the charges^' ^. However, the general form of the field - 
continuous inside the forward light-cone and singular but integrable on the 
forward light-cone - allows to express the retarded propagator by a finite 
number of addends'*: 

GHix,y) = G'^\x,y) + 

because higher orders in the T cannot contribute, as lines of constant X- 
or can only be crossed oncc^^^. There G^^'^{r},y) represents the homo- 
geneous solution in the radiation field only and: 

r^'^(x)+p+'-=p+--(±l)^-.«7' n exp[±za^;«JxT)]. 



-.L,R 



'^The dependence on k is not shown below. 
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V indicates that the above faetors are to be ordered in such a way that, if 
TOi.iJ > ^L^R then > (^-'+)ni„H- The subscript ± stand for the 

sign of the time difference xq — yo and 2p^ = 1 ± 7*^7^. Here, the ultra- 
relativistic Umit has already been taken: In other words, in very energetic 
nuclear collisions the nuclei are highly Lorentz-contracted. For infinitely 
high energies 7 ^ cx), the {b^'^)nR,L go to zero like 7"^. Even now, the 
longitudinal structure is still manifest in the path-ordering of the products. 
However, care has to be taken, as potentially important, because qualita- 
tively different contributions to observables could be neglected, although 
the expression for the propagator in this limit would still be correct. If fi- 
nite longitudinal distances were kept, every charge would have to be treated 
on its own and the previous equation would contain contributions with up 
to Nl + Nr insertions. 

In order to proceed from here, the exact solution for the radiation field 
would be needed which requires transverse lattice calculations beyond the 
scope of this paper. Therefore, the propagator in a model for the radiation 
field A''"'' is to be investigated. By virtue of longitudinally boost-invariant 
boundary conditions on the light cone, it is often taken as being invari- 
ant under longitudinal boosts, too. Although, this is not necessarily the 
case just such a field is chosen here: toA± = ±x±acxp{—x+X-/TQ}, which 
satisfies the Lorenz, Fock-Schwinger, and Aq — gauge-conditions. 

Evaluating Eq.(2) for this field can be interpreted according to I with the 
above gauge of reference. In other words, Wilson lines along the temporal 
direction have to be added to the correlator. The single-particle scattering- 
operators to the lowest order in an expansion in the transverse components 
of the covariant derivatives read Tr = 7+7^ -|- 7-7^ where 



^R ip, Q) = Taro- j- a > ■ j x 



/[r+r-ToV(/x + ^ + 2)] - /b±r^ro7(/i + v + 2)\ 
r±2[(^+l)g^-(i. + l)p^] 

which could be partially resummed into hypergeometric functions. f{a) = 
cr^e~'^Ei(cr) — a and r± = p±+ q±. Note that an expansion in covariant 
derivatives and/or the mass preserves the properties of the propagators 
under gauge transformations. Evaluating the integrand yields here: 

\u{q)Tn{q, -p)v{p)\^ = J2[S{q ■ Tr){p ■ Tr^) - A{p ■ q - m'){TR ■ Tr^)] 

n 

The sum runs over the eigenvalues A„ of the colour matrix a = 
X]„ A„ \n) {n\, where all occurences of a in T are replaced by A„. If a 
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light-cone gauge, e.g. A_ = was chosen as reference, then 7}j = and 



would have to be used. 

The expression for the advanced single-particle scattering-operator 
needed in the case II can be obtained from TA{p,q) = —T^{q,p). For 
the present model, the integration over d*l in Eq.(3) can be carried out 
analytically. However, the result - moreso after carrying out the trace - is 
to lengthy in order to be spelled out here, but will be presented elsewhere^. 
Further, it is not yet clear how the momentum variable k in Eq.(3) can be 
related to the antiparticle momentum p in Eq.(2). More reflection on that 
point is also required. 
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/[r+r-ro V(/x + iy + 2)]- fjp+r^Tq^ /{ji + u + 2)] 
r+^[{li + l)q--{v + l)p_\ 



